Optical nonreciprocity in a piezo-optomechanical system
Xiao Yu-Ming1, Liu Jun-Hao1, Wu Qin3, Yu Ya-Fei1, †, Zhang Zhi-Ming2, ‡
Guangdong Provincial Key Laboratory of Nanophotonic Functional Materials and Devices, and School of Information and Optoelectronic Science and Engineering, South China Normal University, Guangzhou 510006, China
Guangdong Provincial Key Laboratory of Quantum Engineering and Quantum Materials, South China Normal University, Guangzhou 510006, China
School of Biomedical Engineering, Guangdong Medical University, Dongguan 523808, China

 

† Corresponding author. E-mail: yuyafei@m.scnu.edu.cn zhangzhiming@m.scnu.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 61941501, 61775062, 11574092, 61378012, and 91121023), the Doctoral Program of Guangdong Natural Science Foundation, China (Grant No. 2018A030310109), the Doctoral Project of Guangdong Medical University, China (Grant No. B2017019), and the Open Project of the Key Laboratory of Low Dimensional Quantum Structures and Quantum Control of the Ministry of Education, Hunan Normal University, China (Grant No. QSQC1808).

Abstract

We theoretically study the optical nonreciprocity in a piezo-optomechanical microdisk resonator, in which the cavity modes and the mechanical mode are optically pumped and piezoelectrically driven, respectively. For asymmetric optical pumping and different piezoelectrical drivings, our system shows some nonreciprocal optical responses. We find that our system can function as an optical isolator, a nonreciprocal amplifier, or a nonreciprocal phase shifter.

1. Introduction

Nonreciprocal components such as isolators and circulators are indispensable in photonic systems.[1] So far, a number of schemes have been proposed to achieve the nonreciprocity using various physical effects or systems, such as the magneto-optical effect,[2] the interband photonic transitions,[3,4] the stimulated Brillouin scattering,[57] the parity-time-symmetric system,[8] and the rotating nonlinear devices.[9] Furthermore, scientists have found that the nonreciprocity can also be obtained by the optomechanical interaction,[1017] and the realization of optomechanically induced nonreciprocity has been experimentally demonstrated.[18,19]

Cavity optomechanics studies the radiation-pressure interaction between electromagnetic field and the mechanical motion. A large number of physical effects have been demonstrated in cavity optomechanical systems, such as the optomechanically induced transparency,[20,21] the squeezing of light,[22,23] the parity-time-symmetric breaking,[24] the photon blockade,[25] and the cooling of mechanical resonators.[26,27] In recent years, more physical phenomena have been demonstrated in hybrid optomechanical systems, which combine the optomechanical system and other physical systems,[2832] in which the piezo-optomechanical system is attractive since it provides a platform to investigate light-mechanical motion-electricity interaction.[33,34] In the latest studies, the double optomechanically induced transparency,[35] the microwave-optical frequency conversion[36] and the entangled microwave-optical photon pairs[37] in the piezo-optomechanical systems have been demonstrated. In our present work, we study the nonreciprocal transmission, amplification, and phase shift in a piezo-optomechanical system.

This paper is structured as follows: In Section 2, we describe the proposed scheme and give the Hamiltonian of the system and the dynamical equations. In Section 3, we study the nonreciprocal transmission and amplification phenomena. In Section 4, we show the nonreciprocal phase shift. In Section 5, we conclude by summarizing the results.

2. Model and dynamical equations

Our proposed scheme is shown in Fig. 1. We consider a piezo-optomechanical microresonator, in which a clockwise circulating cavity mode a and a counterclockwise circulating cavity mode c couple with the mechanical mode b via radiation pressure. Meanwhile, the mechanical mode is piezoelectrically driven by a microwave field. The mechanical motion of the microdisk is controlled by the piezoelectric effect, which has been described in many works.[3840] The piezoelectric effect in microdisk naturally couples the microwave mode to the mechanical mode by introducing an electromechanical coupling energy given by Hem = ∫ (S · e · E) dV, where S, e and E are the strain filed, piezoelectric coupling tensor and electric field. In addition, cavity mode a is simultaneously driven by a pump field εl with frequency ωl and a probe field εpa with frequency ωpa. Cavity mode c is driven by a probe field εpc with frequency ωpc. The Hamiltonian of the system can be expressed as

Here H0 denotes the free Hamiltonian of the system, with a(c), b and d being the annihilation operators of the clockwise (counterclockwise) circulating cavity mode with frequency ω0, the mechanical mode with frequency ωm and the microwave mode with frequency ωd. H1 describes the optomechanical interaction, with gom denoting the single photon optomechanical coupling strength. H2 describes the piezoelectric interaction between the mechanical mode b and the microwave mode d, with gem denoting the piezoelectric coupling coefficient.[3638] H3 describes the interaction between the cavity mode a, the cavity mode c, the microwave mode d and the driving fields. The amplitudes of the pump field, the probe fields and the microwave pump are defined as , , and , respectively, where Pl, Ppa, Ppc and Pd denote their powers, κe and γe are the extrinsic coupling losses. In Eqs. (1), ϕpa (ϕpc) is the phase difference between the probe field εpa (εpc) and the pump field, and ϕm is the phase difference between the microwave pump and the pump field. The optical mode coupling between a and c results from the imperfections of the disk, for the device we studied, such mode coupling is very small and can be neglected.[38]

Fig. 1. Schematic diagram of our proposed scheme. The piezo-optomechanical microresonator is optically driven by a pump field (amplitude εl, frequency ωl) and two probe fields (amplitude εpa, frequency ωpa, amplitude εpc, frequency ωpc), and piezoelectrically driven by a microwave field with frequency ωd.

In the rotation frame with ħωlaa + ħωlcc + ħωddd, the Hamiltonian of the system is written as

where Δ0 = ω0ωl is the frequency detuning between the cavity mode and the pump field, Δpa = ωpaωl (Δpc = ωpcωl) is the frequency detuning between the probe field εpa εpc) and the pump field. When the microwave pump is strong,[38] the microwave mode d can be treated as a constant number |d|e−iϕm. We take this constant-number assumption, and the effective Hamiltonian of the system can be obtained as

The dynamics of the system can be described by a set of quantum Heisenberg–Langevin equations

where εd = gem|d|, κ = κin + κe is the total loss rates of the cavity modes, with κin denoting the intrinsic loss, γ the loss rates of the mechanical mode, and ain (cin) and bin the δ-correlated operators of the input noises with zero-mean values for the cavity modes and the mechanical mode, respectively.

We have also assumed that εpa, εpc, εdεl, and the microresonator is driven by the pump field at the red sideband (ωl < ω0). Then Eqs. (4) can be linearized by using the ansatz a = α + δa, c = λ + δc, b = β + δb, where α, λ and β are the steady-state mean values, and δa, δc and δb are the fluctuations. By substituting the above ansatz into Eq. (4), and neglecting the nonlinear terms such as δaδb, δcδb ··· and using the rotating-wave approximation,[41] we can obtain the steady-state mean values

and the linearized quantum Heisenberg-Langevin equations

where Δ = Δ0gom(β + β*) = ωm is the effective frequency detuning between the cavity mode and the pump field. We have also assumed that Δpa = Δpc = ωd. By using the ansatz δO = O+e−i ωdt + O ei ωdt (O = a,b,c), we can obtain the following solutions

where δ = Δωd=ωmωd. We can see that the coherent mechanical displacement b+ consists of two terms which correspond to the contributions of the piezoelectric force and the radiation pressure, respectively.

The output fields are given by the input-output relations[42] aout = εl + εpa e−iωdt e−iϕpaκe a, cout = εpce−iωdt e−iϕpcκec. By writing the output field as and , we can obtain

The nonreciprocal transmission is described by the normalized transmission spectra[43]

where η = εd/εpa, ϕ = ϕmφpa, and G = gomα is the effective optomechanical coupling strength. The phase shifts of the output fields are defined as

In the following sections, using above expressions, we numerically study the optical nonreciprocal responses of our system. The parameters are chosen based on a recent experiment:[38] ω0 = 2π × 192 THz, κ = 2π × 1.02 GHz, G = 2π × 0.16 MHz, ωm = 2π × 3.235 GHz, γ = 2π × 0.961 MHz.

3. Piezoelectrically induced nonreciprocal transmission and amplification

From Eqs. (9) we can learn that Ta can be controlled by adjusting the dimensionless amplitude η and the relative phase ϕ. However, Tc is independent of η and ϕ. Now we numerically evaluate Ta and Tc to show the parameter conditions of the piezoelectrically induced nonreciprocal transmission and amplification. Firstly we consider that our system works in the critical-coupled regime (κe = κ/2). As shown in Fig. 2, we plot the transmission spectra Ta and Tc as a function of the detuning δ/κ. When η = 0 [Fig. 2(a)], we find that the system exhibits a reciprocal transmission for the probe fields εpa and εpc. When η is not equal to zero, there will appear a narrow peak in Ta near δ = 0 [Figs. 2(b)2(d)], and the height of the peak increases with increasing η. The height of the peak can be 1 [Fig. 2(c)]. This means a complete transmission of the probe field εpa. In this case our system can act as an optical isolator (Ta = 1, Tc = 0). The height of the peak can even be larger than 1 [Fig. 2(d)]. This means an amplification of the probe field εpa. In this case our system can act as a nonreciprocal optical amplifier (Ta > 1, Tc = 0).

Fig. 2. The transmission spectra Ta and Tc as a function of δ/κ for ϕ = π in the critical-coupled regime (κe = κ/2). The other parameters are (a) η = 0, (b) η = 2, (c) η = 3, (d) η = 3.6.

Then we consider that our system works in the over-coupled regime (κeκ). As shown in Fig. 3(a), the system exhibits a reciprocal transmission when η = 0. When η is not equal to zero, depending on the value of η, Ta will show a narrow valley or peak near δ = 0. Firstly, the valley becomes deeper with increasing η [Figs. 3(b)3(c)]. In Fig. 3(c) we adopt η = 1.5, and we have Ta = 0 and Tc = 1 near δ = 0, i.e., the probe field εpc can be totally transmitted and the probe field εpa is completely absorbed. This means that our system can function as an optical isolator. However, the transmitted direction is opposite with the case described in Fig. 2(c). After this value of η, near δ = 0, Ta will increase with increasing η [Figs. 3(c)3(e)]. When η = 3, the absorption valley of Ta almost vanishes and the system is transparent for both εpa and εpc [Fig. 3(e)]. When η > 3, the system exhibits a nonreciprocal amplification for the probe field (Ta > 1, Tc = 1), as shown in Fig. 3(f).

Fig. 3. The transmission spectra Ta and Tc as a function of δ/κ for ϕ = 0 in the over-coupled regime (κeκ). The other parameters are (a) η = 0, (b) η = 0.3, (c) η = 1.5, (d) η = 2.5, (e) η = 3, (f) η = 3.3.

In fact, the nonreciprocity of our system comes from the asymmetric optical pumping. Because the pump field εl is very strong, which result in the enhancement of the optomechanical coupling between the cavity mode a and the mechanical mode b, and the mechanical mode can be controlled by the mechanical driving, we can indirectly change the transmission of right-moving probe field by adjusting the mechanical driving. On the contrary, the cavity mode c is not driven, and the optomechanical coupling between c and b is too weak to be ignored. Thus the transmission of left-moving probe field cannot be controlled by the mechanical driving and remains unchanged. We note that in a membrane-in-the-middle optomechanical system, the nonreciprocity can also be achieved.[44] Compared with the work in Ref. [44], we show that the nonreciprocal transmission and amplification in a piezo-optomechanical microresonator can be achieved, and we also show that the transmissive direction of the isolator can be changed by adjusting the ratio between κe and κ.

4. Piezoelectrically induced nonreciprocal phase shift

In this section, we demonstrate that our system can perform as a nonreciprocal phase shifter. To achieve a nonreciprocal phase shifter, we should guarantee that both the left-moving probe field and the right-moving probe field are completely transmitted. Figure 4 shows that the probe fields εpa and εpc undergo different phase shifts over a finite range of detuning, and the used parameters are the same as that in Fig. 3(e). This means that our system can act as a nonreciprocal phase shifter. From the application perspective, the nonreciprocal phase shifter can be used to build circulators,[7] which demand the phase difference to be π. Figure 4 shows that |θaθc| can vary from π to 2π, and |θaθc| = π can be met near δ = 0.

Fig. 4. The phase difference |θaθc| as a function of δ/κ. The parameters are the same as those in Fig. 3(e).
5. Conclusion and perspectives

In summary, we have investigated the optical nonreciprocity in a piezo-optomechanical system. The analytic expressions of the transmission spectra and the phase shifts have been given. We have shown the piezoelectrically induced nonreciprocal transmission and amplification in the critical-coupled and over-coupled regimes, respectively. The piezoelectrically induced nonreciprocal phase shift in the over-coupled regime has also been shown. Our works provide a new scheme for achieving the nonreciprocal devices and may have applications in the building of photonic network.

Reference
[1] Bi L Hu J Jiang P Kim D H Dionne G Kimerling L Ross C 2011 Nat. Photon. 5 758
[2] Shoji Y Mizumoto T 2014 Sci. Technol. Adv. Mater. 15 014602
[3] Yu Z Fan S 2009 Nat. Photon. 3 91
[4] Lira H Yu Z Fan S Lipson M 2012 Phys. Rev. Lett. 109 033901
[5] Dong C H Shen Z Zou C L Zhang Y L Fu W Guo G C 2015 Nat. Commun. 6 6193
[6] Kim J Kuzyk M Han K Wang H Bahl G 2015 Nat. Phys. 11 275
[7] Fu W Shu F J Zhang Y L Dong C H Zou C L Guo G C 2015 Opt. Express 23 25118
[8] Arkhipov I I Miranowicz A Di S O Stassi R Savasta S Nori F Özdemir Š K 2019 Phys. Rev. 99 053806
[9] Huang R Miranowicz A Liao J Q Nori F Jing H 2018 Phys. Rev. Lett. 121 153601
[10] Hafezi M Rabl P 2012 Opt. Express 20 7672
[11] Zhang L W Li X L Yang L 2019 Acta Phys. Sin. 68 170701 in Chinese
[12] Xu X W Song L N Zheng Q Wang Z H Li Y 2018 Phys. Rev. 98 063845
[13] Xu X W Li Y 2018 Phys. Rev. 98 053854
[14] Xu X Zhao Y Wang H Jing H 2020 Photon. Res. 8 143
[15] Li B Huang R Xu X Miranowicz A Jing H 2019 Photon. Res. 7 630
[16] Mirza I M Ge W Jing H 2019 Opt. Express 27 25515
[17] Du L Chen Y T Wu J H Li Y 2020 Opt. Express 28 3647
[18] Ruesink F Miri M A Alu A Verhagen E 2016 Nat. Commun. 7 13662
[19] Shen Z Zhang Y L Chen Y Zou C L Xiao Y F Zou X B Sun F W Guo G C Dong C H 2016 Nat. Photon. 10 657
[20] Weis S Rivière R Deléglise S Gavartin E Arcizet O Schliesser A Kippenberg T 2010 Science 330 1520
[21] Yan X Jia W Z Li Y Wu J Li X Mu H 2015 Front. Phys. 10 351
[22] Mancini S Tombesi P 1994 Phys. Rev. 49 4055
[23] Zhang Z C Wang Y P Yu Y F Zhang Z M 2018 Opt. Express 26 11915
[24] X Y Jing H Ma J Y Wu Y 2015 Phys. Rev. Lett. 114 253601
[25] Liao J Q Nori F 2013 Phys. Rev. 88 023853
[26] Liu Y C Xiao Y F Luan X Wong C W 2013 Phys. Rev. Lett. 110 153606
[27] Lai D G Zou F Hou B P Xiao Y F Liao J Q 2018 Phys. Rev. 98 023860
[28] Wang Y P Zhang Z C Yu Y F Zhang Z M 2019 Chin. Phys. 28 014202
[29] Zhang Z C Wang Y P Yu Y F Zhang Z M 2019 Ann. Phys. 531 1800461
[30] Liu J H Zhang Y B Yu Y F Zhang Z M 2019 Front. Phys. 14 12601
[31] Liu J H Zhang Y B Yu Y F Zhang Z M 2017 Opt. Express 25 7592
[32] Liu Z X Xiong H Wu Y 2018 Phys. Rev. 97 013801
[33] Xiong C Fan L Sun X Tang H X 2013 Appl. Phys. Lett. 102 021110
[34] Balram K C Davančo M I Song J D Srinivasan K 2016 Nat. Photon. 10 346
[35] Wu S C Qin L G Jing J Yan T M Lu J Wang Z Y 2018 Phys. Rev. 98 013807
[36] Zou C L Han X Jiang L Tang H X 2016 Phys. Rev. 94 013812
[37] Zhong C Wang Z Zou C Zhang M Han X Fu W Xu M Shankar S Devoret M H Tang H X Jiang L 2020 Phys. Rev. Lett. 124 010511
[38] Fong K Y Fan L Jiang L Han X Tang H X 2014 Phys. Rev. 90 051801
[39] Han X Fong K Y Tang H X 2015 Appl. Phys. Lett. 106 161108
[40] Shen Z Han X Zou C L Tang H X 2017 Rev. Sci. Instrum. 88 123709
[41] Aspelmeyer M Kippenberg T Marquardt F 2014 Rev. Mod. Phys. 86 1391
[42] Gardiner C Collett M 1985 Phys. Rev. 31 3761
[43] Liu J H Yu Y F Zhang Z M 2019 Opt. Express 27 15382
[44] Li Y Huang Y Y Zhang X Z Tian L 2017 Opt. Express 25 18907